Dipolar spin correlations in classical pyrochlore magnets 
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We study spin correlations for the highly frustrated classical pyrochlore lattice antiferromagnets 
with O(N) symmetry in the limit T — » 0. We conjecture that a local constraint obeyed by the 
extensively degenerate ground states dictates a dipolar form for the asymptotic spin correlations, 
at all TV 7^ 2 for which the system is paramagnetic down to T = 0. We verify this conjecture 
in the cases N — 1 and TV = 3 by simulations and to all orders in the l/N expansion about the 
solvable N — oo limit. Remarkably, the N — oo formulae are an excellent fit, at all distances, to 
the correlators at TV = 3 and even at TV = 1. Thus we obtain a simple analytical expression also for 
the correlations of the equivalent models of spin ice and cubic water ice, In- 

PACS numbers: 



Introduction: In frustrated magnets, magnetic order is 
strongly suppressed compared to expectations from sim- 
ple mean-field theory. This gives rise to a cooperative 
paramagnetic, or spin liquid, regime, where the energy 
scale of interactions exceeds that set by the temperature, 
and nonetheless no long-range order ensues; if frustration 
is particularly severe, magnetic ordering can be absent 
altogether. A particular case in point is the pyrochlore 
lattice, where the absence of ordering for the Ising an- 
tiferromagnet was already argued in 1956 1J. The Ising 
model is of particular interest on account of its T = 
entropy and its equivalence to cubic ice Ih, and it has 
also been approximately realised experimentally in the 
titanate spin ice compounds 0- There is strong evidence 
that the Heisenberg model on the pyrochlore lattice does 
not order down to T = 0, 01- an( A a general the- 
ory based on constraint counting indicates that this is 
so for all O(N) magnets with N > 3 although not for 
N = 2, where thermal fluctuations lead to collinear or- 
dering ( "order by disorder" ) [ 




FIG. 1: The pyrochlore lattice. 

The absence of long range order does not mean the 
physics is trivial at low temperatures as the accessible 
states have non-trivial local constraints on them — hence 
cooperative paramagnet. For the thermodynamics it is 
possible to make do with small clusters, e.g. Pauling's 
estimate of the entropy of ice works rather well and the 
thermodynamics of the Heisenberg pyrochlore magnet is 
well-described even by an approximation based on an iso- 
lated tetrahedronQ. For the spin correlations however, 



the situation is presumably different. Indeed, Monte 
Carlo simulations of the Heisenberg magnet by Zinkin et 
aL[8| have demonstrated the presence of sharp features in 
the structure factor in certain high symmetry directions 
(the "bow-ties" in Fig.|21), indicative of extended spin cor- 
relations. It is understood, qualitatively, that these are 
a direct consequence of the ground state constraint that 
each tetrahedron separately has vanishing total spin[{|. 

In this paper we provide a theory of the T — ► spin 
correlations of nearest neighbor O(N) magnets on the 
pyrochlore lattice when they are paramagnetic down to 
T = 0. This restriction excludes only the case N = 2, but 
even there our results are probably applicable as we dis- 
cuss below. Our work builds on two previous advances. 
The first is the recognition RjllfJ, [ll], u3> generalizing 
work on two-dimensional ice|13|. that d = 3 models with 
binary degrees of freedom and a local "Gauss's law" con- 
straint should generally exhibit asymptotic dipolar corre- 
lations governed by a pure Maxwell action. This applies 
directly to our N — 1 (Ising) case, while the underlying 
argument has been conjectured to apply also to Heisen- 
berg spins by Henley [lOj - The second is the solution by 
Garanin and Canals of the N = oo limit , in particu- 
lar their numerical determination of the structure factor 
in the N — oo limit, which exhibits the same bow-ties 
observed by Zinkin et at 

In the following we (a) motivate how dipolar cor- 
relations arise in pyrochlore magnets from local con- 
straints, (b) demonstrate their existence analytically in 
the N = oo limit and their persistence to all orders in 
the l/N expansion, and (c) verify by simulation that the 
correlations at N = 3 and N = 1 are dipolar. We find 
that the correlations at N = 3 are very well fitted by the 
N = oo formulae even without modifications at short dis- 
tances. Remarkably we find that the N = oo formulae 
are probably even more accurate for N = 1 where the rel- 
ative and absolute error never exceed 2% and 4.5 x 10~ 4 , 
respectively even though this is not a conclusion that 
one would guess on the basis of the l/N expansion! We 
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FIG. 2: Top: The [hhk] structure factor for N = oo [dark=low 
intensity; the x (y) axis corresponds to h (k)]. Middle: Ising 
correlations at T = (bottom: Heisenberg, at T/J — 0.005) 
from Monte Carlo simulations, compared to T = large- 
N correlations: (Si(x)Si(O)), for two inequivalent directions, 
[101] and [211], multiplied by the cube of the distance, x 3 . 



The error bars are around 5 x 10" 
curves are for system sizes L = 
There is no fitting parameter. 



(5 x 10' 
8,16,32,48 (L = 8,16,32). 



close by noting the applicability of the ideas in the paper 
to related syst ems, which will be detailed in a separate 
publication 15], as well as their known limitations. 
Local constraints and dipolar correlations: We 
consider the classical nearest neighbor antiferromagnet 
on the pyrochlore lattice, H = j\ SjSj-, where 

the Si are TV-component spins of fixed length N. We 
set J = 1 in the following. The pyrochlore lattice 
consists of corner sharing tetrahedra, and the Hamil- 
tonian can be rewritten, up to a constant, as H = 
(1/2) 2s(SieH Si) 2 ; where the sum in parentheses runs 
over all four spins at the corners of a given tetrahedron, 



B, and the outer sum is over all tetrahedra. Hence, the 
ground states (minimum energy configurations) are such 
that for each tetrahedron and each spin component a, 



S? = 0. 



(1) 



This can be turned into a manifest conservation law on 
the dual - bipartite diamond - lattice, the sites of which 
are given by the centres of the tetrahedra while the spins 
sit at the midpoints of its bonds. 




FIG. 3: The centres of the tetrahedra define (the two sublat- 
tices of) the diamond lattice, denoted by black and grey dots. 
The orientation of the diamond bonds, and the four sublat- 
tices of the pyrochlore lattice, are defined by the vectors e K 



First, we orient each bond, k, by defining a unit vector 
e K , which points along the bond from one sublattice to 
the other, see Fig. [3] Next we define N vector fields 
on each bond, B" = S%e K , where S% denotes the spin on 
bond k. The ground state constraint (y implies that each 
B a separately forms a set of solenoidal fields at T = 0, 
V • B a = 0. 

For N = 1, spin flips connecting two ground states 
correspond to reversing the direction of a closed loop of 
"magnetic flux", B; evidently, B averages to zero over 
such a flippable cluster of spins. Upon coarse-graining, 
a high density of flippable clusters (and therefore a large 
number of ground states) translates into a small (well- 
averaged) coarse-grained B. We now posit that this fea- 
ture carries through to N > 1, so that states with small 
values of B a will in general be (entropically) favoured. 
This is captured most simply by introducing a weight 
functional, p: 



P [{B Q (x)}] ex exp 



K 

y 



(2) 



provided the solenoidal constraint is implemented, where 
K is the stiffness constant. (If we solve the constaint by 
introducing a vector potential for each component, we 
are led to the Maxwell action.) 

From this we can deduce the long distance correlators, 

r 2 Sa 



(Sf (x)Bf (0)) cx 6*0 



3 X j tjj n 



(3) 
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which are dipolar as advertised. The appropriate long- 
wavelenth formulae for the case N = 1 have already been 
given m Ref.H 

This argument does not take into account thermal fluc- 
tuations out of the ground state manifold. These are 
gapped, and thus unimportant for T — ► 0, for N = 1. 
However, for continuous spins they endow each micro- 
scopic ground state with a non-trivial entropic weight, 
and are essential for even defining what is meant by a 
measure on the set of ground states. An analysis based on 
Maxwellian constraint counting implies that this weight- 
ing is sufficiently uniform to be innocuous for N > 3 
so that the dipolar forms should hold. For N — 2 this 
is not the case and the weighting leads to an order-by- 
disorder phase transition 0, • However at this tran- 
sition, the spins order collinearly but do not pick fixed 
orientations along their common axis. Hence, for the 
resulting collinear ensemble, the N = 1 results can be 
applied, with the caveat that it has not been established 
whether a further phase transition will select a subset of 
collinear states at even lower temperatures. 

While the arguments given above are based on the 
plausible ansatz (0) , we now turn to establishing their ac- 
tual correctness. The continuity of the physics for N > 3 
will allow us to tackle this in the 1/N expansion and we 
will supplement this with explicit simulation at small N. 
N = oo: We start with the simple classical O(N) model; 
its N = oo saddle point^(| is, of course, the one dis- 
cussed in Refs. 14. The partition function is given by 
Z = f T>cf)D\exp[— S(tfi, A)], with the action defined as 



N 



i } j a—1 



where Ay is the interaction matrix divided by T. 

The saddle point evaluation for N — oo is analogous 
to that of Refs. as we are n °t aware of an analyt- 
ical treatment of the full pyrochlore correlations in the 
literature, we provide it here. 

The Fourier transform of the interaction matrix is 
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where c (, = cos( 

values of the interaction matrix are |17| 

fi,2 = 0, V3A = 2 ± y/l + Q, 



->lb 



(4) 



). The eigen- 



where Q = ^ + (^ + (4 + 4^+^ + ^-3. 

These can be combined to give the spin correlators [f4j. 
For the zero temperature case discussed here, only the 
modes with v = contribute. There are two independent 
correlators, namely those between spins on the same and 
between spins on different sublattices. We thus find the 



following spin correlators (we label four sublattices as 
shown in Fig. |3J 



6 - 2r 2 - 2c 2 - - 2c 2 - 
(S?(q)Sf(-q)) = 5 afi xz yz 



3-Q 



(S?(q)S£(-q)) = 



2(cos(f) 



2 J^xz *~-xz } 



3- 



(5) 



From these one can calculate the correlators of B" = 
&k c« an d verify that they have the small q forms, 

q 2 



(Bf(q)^(q))«M««-^) 



(6) 



equivalent to the real space forms |J5} noted earlier. 
The structure factor in this limit is 



S(q) = 4 
+ (C 



(c yz + Cyz) sin 2 (Jj) 



3-Q 

xz + c—) sin 2 + (c xy + c^y) sin 2 (J^j 



•(7) 



It has been customary to consider the structure factor 
of the pyrochlore antiferromagnets in the [hhk] plane in 
reciprocal space, as it contains the high-symmetry di- 
rections [100], [110] and [111]. Here, the wavevector 
q = and Eq.0 simplifies to 



3 2(cos(f)-cos(f)) 2 sin(f) 2 
5 — cos(q x ) — 4cos(^) cos(^) 



Si 



[hhk] 



(8) 



This is plotted in Fig. |3 the same plot has previously 
appeared in Ref . Q. We emphasize that it is the trans- 
verse, dipolar form of the correlators which gives rise to 
the bow tie structures. 

1/N expansion: To set up the 1/N expansion [l^, we 
expand away from the saddle point by allowing the La- 
grange multiplier field A to vary from its uniform saddle 
point value, A: Aj = A + /i^. Integrating out the <j>'s 
removes the term linear in /i; the quadratic term yields 
the propagator and higher order terms generate vertices. 
The quadratic term, J2ij Vii&^j) 2 Mi = ■ ViM^ 1 Hj , 
implies a propagator (/ij/ij) = Aiij- One can argue, 
and we have checked by explicit numerical computation, 
that (l^i^j) decays inversely as the sixth power of the 
separation between i and j. Hence its Fourier trans- 
forms, M KK '(q), are continuous functions over the Bril- 
louin zone. We note that the N = oo correlators of the 
spins G Q KK , (q) are discontinuous at q = on account of 
their dipolar character. 

We now briefly sketch the proof that the dipolar form 
on the spin correlations survives to all orders in the 1 /N 
expansion; details will be presented in . Consider the 
Dyson series for the spin correlator, G, at finite N: 

GV(q) = GL'(q) + G° KKl (q)S KlK2 (q)G° 2K , (q) + ■ • ■ 

where the self energy E is the sum of the relevant graphs 
at all orders in 1/N. For each graph that enters this sum 
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one can show (I) that £fcz(q) is a continuous function of 
q in the Brillouin zone, (II) that lim q ^ (q) = a dlag 
for k — k' and a od otherwise. (I) follows from the fact 
that in all such diagrams one integrates over an integrand 
that is bounded and at most discontinuous at isolated 
points, while (II) uses the action of lattice symmetries on 
G° and M. . Using these observations the leading small q 
terms in the Dyson series lead to a simple geometric sum 
and one finds that 

GW(q)cxG^,(q) (9) 

near q = 0. Consequently, the correlations are still dipo- 
lar at long distances. However, even in the most opti- 
mistic case we would not expect the series to converge 
beyond N = 3, where long range order sets in due to 
order by disorder. 

Comparison to Monte Carlo simulations: To the 

naked eye, the N — oo structure factor of Fig. [21 looks 
almost identical to the one obtained by Zinkin for Heisen- 
berg spins0 and by us for Ising snins[lij. For a more 
sensitive test we have compared the real-space expres- 
sions to Monte Carlo simulations for Ising and Heisenberg 
spins for different system sizes, shown in the same figure. 
We find that the error for Ising spins is very small indeed 
- the relative and absolute systematic error never seems 
to exceed 2% and 4.5 x 1CT 4 respectively. The agreement 
for Heisenberg spins is somewhat worse. However, since 
the disagreement is of the size of the finite-temperature 
error bars, we cannot extrapolate down to T — > reli- 
ably. In either case, one sees that the N — oo formu- 
lae we have derived work remarkably well quantitatively 
even without any further corrections — indeed, far better 
than one might have guessed at the outset — even at small 
distances, and they correctly capture finite-size effects. 
Applicability to other models: As the Ising anti- 
ferromagnet on the pyrochlore lattice is isomorphic to 
(nearest- neighbour) [l^ spin ice, the results can straight- 
forwardly be transcribed, and further carried over to cu- 
bic (water) ice 7/j, which is again equivalent 0. In the 
latter context the asymptotic form of the correlations is 
known 0] and the interest of our results is in the accuracy 
of our analytic forms at all distances. The analysis pre- 
sented here straightforwardly generalises to other three 
dimensional cooperative paramagnets with conservation 
laws. 

The present approach can also be applied to two- 
dimensional models with a conservation law. However, 
it breaks down for models with a discrete set of ground 
states. Technically, the effective height (two-dimensional 
gauge) field acquires a non-zero compactification radius 
in this case, which leads to the appearance of additional 
operators under coarse-graining. These are absent for 
N = 00, and their appearance is non-perturbative in 
N. For the kagome magnet, for example, they arise for 



N = 2,3, which is why the y/3 x V3 correlations known 
to be present there were not found in Ref. 14. However, 
information gleaned from the present approach, such as a 
stiffness, can nonetheless usefully be fed into these mod- 
els. Finally, we have shown that the present approach 
can also be applied to a class of paramagnetic models 
without a conservation law, for which so far very little 
has been known in terms of correlation functions. This, 
along with the more technical material, will be discussed 
elsewhere I la. 
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